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I. IN IBCEUCTION 



A variety of acoustic models exist to predict transmis- 
sion Jess. Each of these models contain inherent strengths 
and weaknesses. All have shown poor results in a shallow 
water environment due to difficulties at the ocean fcettem 
where there is an interface tetween media of different sound 
speeds and densities. 

Since its introduction (Hardin and Tappert, 1 9 7 3 j , the 
parabolic wave eguaticn has been a widely accepted means cf 
solution for acoustic propagation. The earliest programs 
used a split-step Fourier transform algorithm to solve tne 
parabolic eguation (PF) . Several other solution techniques 
have beer developed primarily to overcome difficulties that 
occur when the Fourier transform encounters an interface 
between different media (Lee and Botseas, 1982 and McDaniel 
and lee , 1 9 82) . 

An alternative solution technigue that uses an implicit 
finite-difference (III) algorithm was developed by Lee and 
fapacakis (1979). Ibis method incorporates appropriate 
interface conditions and allows solutions in shallow water. 
Starting with the IFE algorithm, Jaeger (1983) developed a 
computer mcdel to predict transmission loss and acoustic 
pressure based on user specified bottom topograph} ar.d a 
single sound speed profile. This computer model uses the 
mathematical treatment of the horizontal and sloping inter- 
faces developed by McDaniel and Lee (1982) and lee and 
McDaniel (1983). It also utilizes several design features 
and methods incorporated in an earlier computer program 
developed by Lee and Eoteas (1982), and a PE computer model 
developed by 3rcck (1578). 
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The ifD program preserves continuity of pressure and 
continuity cf the ncimal component of particle velocity at 
an interface between media having different sound speeds and 
densities. This feature makes the program unconditionally 
stable and better able to handle the bottom boundary 
condition. 

Since its development, the IFD program has net been 
rigorously tested. This thesis analyzes the program's 
perf crmarce in an idealized shallow water environment. The 
environment includes a simple sloping sand bottom beneath ar. 
isospeed water field. The analysis begins by comparing the 
IFD ' s predictions with predictions from Jensen and 
Kupernan's { 1 9 8 0) El model and Coppens, Humphries and 
Sander's (1S84) image model. The analysis also includes a 
comparison of the ncael' s estimated transmission less 
contours with expectations based on simple physical 
reasoning. Finally, the thesis describes an attempt to model 
a shallow sloping bottom in an experimental tank. The tank 
contains a sand bottom sculptured with a ten degree slope, 
laboratory measurements cf tne pressure field are taken at a 
frequency cf 100 kH 2 for comparison to the predicted pres- 
sure field generated by the IFD. These comparisons cf the 
IFD predictions with ether model estimates, theory, and 
laboratory measurements give an indication of the IFE's 
performance in a shallow water environment. 



II - THE IFE CCMPOTER MO DEI 



1. EACKC-PCOND 



In implicit finite-difference solution technique tc the 
parabolic equation has been studied and refined by many 
authors. The history of this development is explained in 
detail by Jaeger {1382 } , but merits review in order tc gain 
a perspective on an analysis of the IFD computer model. The 
parabolic eguation is an approximation to the elliptical 
wave eguaticn. The first means of solving the PE used a 
split-step fast Fourier transform method as developed by 
lappert and Hardin (1973). This method requires periodic 
boundary conditions ir depth because of the finite Fourier 
transform and handles this constraint by introducing an 
artificial horizontal pressure release bottom belcw tne 
actual physical bottcm. This method of implementing an 
artificial bottom was incorporated into the earliest PE 
models developed by Jenson and Krol ( 1975) and Brock (1S76). 

Errors in this split-step Fourier transform method were 
found tc be proportional tc the horizontal range step and 
the second derivative of the index of refraction. The second 
derivative of the index of refraction tends to be large 
across the ocean bottcm interface. Another problem with the 
split-step method is that it does not consider density 
differences between two different media at an interface, 
which influences the reflection coefficient. For these 
reasons the split-step fourier transform method proved to be 
poorly suited for a shallow water environment. 

The IFE soluticr method was introduced by lee and 
Papadakis (1979) as ar alternative to the split-step method. 
The IFD method employs a second order central difference 
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formula to solve the EE in the form ox a tridiagonal matrix, 
although the first version of the IF D did handle discontinu- 
ities in the sound speed profile, it did not consider 
density discontinuities. In 19 62 idcDaniel and Lee introduced 
a method to handle a horizontal interface of different 
densities. In 1983, they extended tneir treatment to 
include a sloping interface. It is this version of the IFD 
that is used in Jaeger's computer program. 

£. IEI CChEOTEE HOD II. 

Ihe 1FE computer program consists of a main program and 
twenty subroutines. Ihe pregam utilizes a modular construc- 
tion sc that each of the various subroutines are called from 
the main program to complete a specific calculation or func- 
tion when required. Ihe IFE is run interactively from a 
user generated input file that contains values for 
freguency, one sound speed profile, a bottom profile, 

source/receiver depths, attenuation coefficients for betn 
the water and the bottom, and several other input parameters 
that tell the the program where to obtain a solution within 
the field. The program initiates the calculations assuming 
an initial Gaussian pressure field and an artificial pres- 
sure release surface at a user specified depth. 

Attenuation in both the water and tne sediment is 
handled as complex indices of refraction. An artificial 
attenuation layer is established beneath the sediment to 
introduce attenuation above the artificial pressure release 
surface. The actual magnitude ox this enhanced attenuation 
is calculated using ax eguation derived by Brock (1978) for 
use in his EE computer model. 

The IFD program steps along the specified bottom profile 
and makes calculations down tnrough the wat ex/sediment 
column at each user specified horizonal range. The program 

/ 
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C. eCEEI IFOBLEHS/MCIIFICAIICNS 

As this study of the IFD program progressed, it became 
necessary to modify certain aspects of the model. Most of 
these modifications were necessary to alter the program 
output irtc a mere desirable form, but a few were imple- 
mented to correct picgramming deficiencies. Although this 
section cf the thesis discusses the earliest model runs, 
these results are ret presented in detail, but orly 
discussed in general terms because they were obtained before 
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acoustic field. There was no observable decrease in trans- 
mission less at the various dump distances as expected. 
Bather, results indicated widely fluctuating patterns in the 
acoustic rield that appeared inconsistent with both previous 
studies and simple physical reasoning. Upon closer analysis, 
it was discovered that although tne program was designed to 
he independent of scale there are several logic statements 
that are net implemented if the user provided range step is 
less thar one meter. Because the logic statements aren't 
satisfied, the h EWSZG and N EjiMAT subroutines (Jaeger, 1SS3) 
are net called correctly. The NEKSEG subroutine initialises 
the bottom slope anc the NEWMAT subroutine computes matrix 
elements for the program. Obviously, errors in these two 
program functions seriously distort results. Because of this 
systematic error in the program it became necessary to scale 
up all tank parameters. All distances were scaled up by a 
factor of 1C00, and frequency was scaled down by a factor of 
1000. Careful analysis reveals that all input parameters are 
a function cf either distance or frequency, so this scaling 
produces results that model those expected in the tank. 

The second modification cf the If D was required- to 
provide a three dimensional graphics display. As discussed 
earlier, the IFE provides a transmission loss plot versus 
range at a single depth. However, to study the model 
predictions in greater detail it was felt that a two dimen- 
sional analysis cf the model estimates would be mere mean- 
ingful. As a result a transmission loss contouring program 
was developed. The program (Appendix B) displays transmis- 
sion less contours for range versus depth. Use cf the 
contour plot requires that transmission loss values gener- 
ated by tie IFD be sent to a data disk used by the contour 
routine. To facilitate this transfer a dummy variable (EltD) 
was established in the PRINT2 subroutine to store the trans- 
mission loss values and then these values are written to the 
data disk at the end cf the main program. 
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Another modification of the IFD output was required to 
change the real and imaginary components of the pressure to 
a single pressure amplitude magnitude. It was felt that 
dealing with the pressure magnitude was easier ard mere 
meaningful than with the real and imaginary components of 
the pressure field. lhis conversion was done in the PEINT2 
subroutine and established a new variable (PRMAG) to repre- 
sent the pressure magnitude. 

A final modification in the computer program was made 
due tc a suspected error in the computation of the attenua- 
tion in the artificial layer. Physical reasoning dictates 
that proper implementation of the artificial attentuation 
would result in a steady drop off in acoustic pressure with 
depth throughout the artificial layer# with pressure drop- 
ping tc reio at the pressure release surface. IfC model 
results on the other hand shoved wide fluctuations in pres- 
sure vith depth in the layer and then only a minimal fall 
off at the pressure release surface. The equation in the 
IFD that actually computes the magnitude of the attenuation 
in the artificial layer was taken directly from Ercck's 
< 1 9 7 € J PI model (Jaecer, 1S83). However, Brock’s equation 
was derived with feet as the unit of measurement while 
Jaeger's model is derived with meters as the unit of meas- 
urement. With this in mind, Jaeger's equation should be 
approximately a factor of three larger than Brock's eguation 
to correct for the difference in units. To correct for this 
error the equation tc calculate attenuation (ATT (I)) in the 
NEWMA1 subroutine was increased by a factor of three. When 
this correction was implemented, the large fluctuations in 
pressure with depth were eliminated. The expected drop off 
in pressure with depth and the fall off of pressure tc zero 
at the pressure release surface were noted. 

A listing of the revised IFD computer program with all 
modifications can be seen in Appendix A. 
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HCE El V ERIF I CATIC F 



1 • Cc c p ari sen with Jaeger Model Ba r 

In light of the modifications to the I FD computer 
program just discussed, it was necessary to ensure that the 
the dances themselves did not introduce errors ntc the 
model. So as a first step the modified IFD program was run 
for one of the test cases used by Jaeger in his original 
work. His case analyzes propagation in an environment that 
moves from deep to shallow water. This particular environ- 
ment is depicted in Figure 2.1 and was chosen because it was 
very similar to the simple sloping bottom in the tank exper- 
iment. A solution is obtained for a bottom with an upsiepe 
of 8.5 decrees. An isospeed water field is used with sound 
speed set at 1500 m/s. Source frequency is 25 Hz. Ihe 
scenario has a maximum depth of 350 meters and a range of 40 
kilometers. Both the source and the receiver are set at 25 
meters. 

Ihe results using the modified IFE program were the 
same as Jaeger's for this test case above the artificial 
attenuation layer. As explained earlier, the program was 
modified to reflect higher attenuation in the artificial 
layer, so as to properly reflect the effect of attenuation 
in this region. The modified program results did net show 
the large fluctuations in pressure with depth in the artifi- 
cial layer but ratter the gradual decline in pressure 
towards a value of zeio at the pressure release surface. Eut 
above the artificial layer the modified program results were 
identical with those achieved by Jaeger with the original 
program. Apparently, the miner changes in the program 
designed to improve on the form of the program output dees 
not hinder the model's ability to achieve a solution in the 
upper sediment and water column. 
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figure 2. 1 Jaecer's Deep-to- Shallow Water Case. 

2. Cc m c ar i s cd with Jensen and Ku p er man M odel Run 

lie second attempt at verifing the IFE model 
involved a comparison of model results with those achieved 
with a ff model designed by Jensen and Kuperman (1980). Ihe 
Jensen and Kuperman Model (JKM) uses a split-step solution 
technigue. Comparison with this particular model was chosen 
because it is one of the few that obtains a solution in two 
dimensions. Most acoustic models obtain a solution at only a 
single depth. In addition, Jensen and Kuperman made their 
model runs in a simple sloping ocean bottom environment very 
similar to the scenario of interest modeled in the tank. 
Ihis environment is depicted in Figure 2.2 and features a 
gently sloping bottom of 2.2 degrees. Ihe water column has a 
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uniform speed of 150C m/s. The source is placed just helow 
the midpoint in the channel at 112 meters and is driven at a 
frequency cf 25 Hz. The maximum depth in this scenario is 
200 meters and has a maximum range of 12.5 kilometers. 

The Jensen and Kuperman results were expressed as 




.Figure 2.2 Jenser and Kuperman Sloping Bottom Case. 

transmission loss contours fcr range versus depth. Their 
study concentrated cn transmission loss patterns in the 
sediment, hut results were obtained both in the sediment and 
in the water. These results are compared to the IFD results 
at ranges cf 2.5, 5.C, 7.5, 10.0, and 12.5 kilometers. These 
ranges were chosen fcr analysis because the JKM results 
showed the greatest variation in transmission loss with 
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depth and thus, makes for a more meaningful comparison with 
the HE results. 

lie IFD and JKb results can be seen in Figure 2.3 
through Figure 2.7. In all the figures, the IFD estimates 
are shcwr as a solid curve while the Jensen and Kuperman 
results are depicted as circular points. Ihe first analysis 
is at a ranee of 2.5 kilometers. From Figure 2.2, it can be 
seen that the depth at this range is 200 meters and is ir a 
flat lettem region. From the results shown in Figure 2.3, 
it is ctvicus that loth models obtained almost identical 
results from the sirface down to a depth of about 300 
meters. Ihere are differences between the two sets or 
predictions from the ocean bottom to a depth of 100 meters 
lelow this point. Eelow 300 meters the Jensen and Kuperman 
results shew a very slight increase in transmission less 
(IL) with depth. The IFD also shows an overall increase in 
transmission loss with depth but with several fluctuations 
in transmission loss and a marked peak at about 300 meters. 
So in general, the results from the two models have the same 
general tendencies although the IFD appears to shew greater 
detail in results near the water/sediment boundary. 

Figure 2.4 shews results at 5.0 km in range. Here, 
the water depth is still 300 meters and marks the very 
beginrirg of the sloping bottom section. For this range the 
JKM predictions are only available to a depth of 300 meters. 
Ihe results are nearly identical with those obtained by the 
IFD. Both models show a relative minimum in transmission 
loss at a depth of ICO meters and then a gentle increase in 
11 with depth. 

The model results at a range of 7.5 km are seen in 
Figure 2.5. At this range the bottom depth is aheut 150 
meters and the bottom is sloping. Ine models shew the 
greatest difference at this range. Both models produce 
nearly identical results in the water column, but beneath 
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Figure 2-3 IFC and JKM Comparison at a Bange of 2-5 Km. 
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Fignre 2.4 IFD and JKK Conjcarison at a Bange of 5.0 Km. 
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Figuze 2.5 IFE and JKM Comparison at a Bange of 7.5 Km. 




Figure 2.6 IFE and JKH Comparison at a Range of 1C.0 Ka. 
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figure 2.7 IPE and JKfl Comparison at a Bange of 12.5 Km. 

2 6 



150 meters in depth the model estimates iegi n to shew 
differences. The Jensen and Kuperman results show a gradual 
increase in TL with depth below the ocean bottom. The 1FD 
results cn the other hand, shew a gradual decrease in trans- 
mission less from 15C meters to 300 meters, with a relative 
minimum at 300 meters and then an increase in TL beneath 
this depth. Although there are some differences in model 
estimates near the boundary in the flat nottom sections 
previously discussed, the differences appear to be greater 
at the boundary in this sloping bottom case. It these 
regions away from the boundary in either the water column or 
in the sediment, however, both models produce similar 
results. The JKfl's difficulty in obtaining an accurate solu- 
tion at the sediment boundary is not totally unexpected. Ine 
model uses a variation of the straight split-step solution 
technigue similar tc Brock's computer model (Jenser and 
Kuperman, 1S80) , that has characteristically been unsuc- 
cessful in obtaining a reliable solution near a boundary. 

The results at a range of 10-0 km are seen in Figure 
2.6. At this range the water depth is approximately 100 
meters and the bottom is again in a sloping region. At this 
range Jensen and Kuperman results are only available to 
about 15C meters in depth. For the data available, the 
models produce nearly identical results. Beth models predict 
a minimum at about 5C meters in depth and then an almost 
linear increase in 11 with depth. The IBB results also 
reflect c transmissicr loss maximum at about 200 meters, and 
then a slight decrease of II beneath this depth. Bata at 
these depths are not available from the Jensen and Kuperman 
run. Ir light of tie results obtained at 7.5 km for a 
sloping bottom case one might expect that the results for 
the two models would show differences near the ocean bottom. 
However, since there is cnly one Jensen and Kuperman 
prediction available rear the bottom for this range, it is 
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difficult tc make any definite conciusicns regarding differ- 
ences ir model performance at the boundary. 

Eesults at tie apex (range equal to 12.5 k re ) can be 
seen in figure 2.7. Since all estimates are made in the 
sediment at this range, there is no boundary tc contend 
kith. At this range the two models snow the test agreement. 
Eoth models show a gradual decrease in transmission less to 
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influenced by a boundary. In both the water and the deep 
sediment the two models produce similar results. Although 
this is ret conclusive evidence, these similarities suggest 
that the IPX can successfully model acoustic propagation in 
these regions. Near the water/sediment boundary however, the 
Jensen and Kuperman aid If D predictions show marked differ- 
ences. Ihese differences appear to intensify as the hettom 
becomes mere sophisticated. In general, the Jensen and 
Kuperman results do ret seem to show the detail the IPD 
results do. Considering the different solution techniques 
employed by the two models, these differences in results are 
expected. 

2. Com p ari sen with Cop pens, H umphri es and Sanders ho del 
Pun 



Ihe third attempt at verifying the IFD performance 
was dene by comparing results with an image theory modei 
derived by Coppens, Humphries, and Sanders (1984). This 
Coppens, Humphries, and Sanders Model (CHSM) uses a saddle 
point approximation to an image model tc solve for the 
acoustic field. Beth programs were run for the scaled 
scenario modeled in the tank experiment. This scenario is 
depicted in Figure 2.8 and features a ten degree sloping 
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bottom. Ihe maximum water depth for the run is 350 meters 
and the maximum range is two kilometers. The source gener- 
ates a 100 Hz signal and is set at a depth of 175 meters. 
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Figure 2.9. In the Figure# the 
the solid curve while the image 
circular points. Ihe predictions 
malized pressure amplitude for a 
models produce pressure amplitude 
and thus# had to he normalized to 
possible. This normalization was 
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cone t cr each model ly dividing the presssure amplitude by 
the maximum amplitude predicted by the model. 

Ihe results suggest that even though the models 
predict similar larce scale trends in pressure amplitude 
with depth, there are several differences in detail. The 
plot shews that both the IFE and the image model predict an 
almost linear increase in amplitude with depth down tc a 
specific maximum and then a slow decrease in amplitude 
beneath this maximum. The HE however, shows tne maximum at 
about 15.5 meters in depth while the image model maximum is 
deeper at IS meters. Beneath this maximum the IFD shews a 
much sharper decrease in amplitude than the image model. The 
trend in the image mccel data appears smooth, wnile the IFD 
curve reflects several small scale fluctuations. 

It appears that both models predict similar large 
scale trends in pressure amplitude for this scenario. It is 
difficult, if not impossible, to account for the differences 
in detail. As a minimum, however, this comparison indicates 
that the two sets cf predictions are consistent with cne 
another and can be considered reasonable in this shallow 
water environment. 

h . Comparison with Phy sical Seasoni ng 

Frca model comparisons it appears that the IFE at 
least makes a reasonable estimation of acoustic fields ir a 
shallcw water envirorment. IFD results are also analyzed in 
comparison with basic physical reasoning and theory as a 
fourth attempt at model verification. This analysis examines 
IFD tr ansmi ssion loss contours for the simple ten degree 
sloping ccean scenario seen in Figure 2.8. The scenario is 
two kilometers ir ranee with a maximum depth of 350 meters. 
The the source generates a signal at 100 Hz and is placed at 
175 meters in depth. 
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ligore 2.9 IF! and CESfl Comparison at the Apex. 
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SIN ( khd/r) = 0 



cr in ether words; 
khd/r = n^ 



where: 

n 0,1, z . . . 

r = Hange 

X = Wavenumber (2^*/,)) 
d = Scurce Depth 
h = Depth of Node 
^ = Wavelength. 
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Figure 2,10 IFD Tt Contours (dB) from the Source to 604 Meters 
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Figure S.11 IFD TL Con to urs (db) from 600 tc 1350 Meters 
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Figure 2.12 IF D TL Contours (db) froi 1350 to 2100 Meters 



o 







36 



Figure 2.13 IFD TL Contours (db) fro* 1550 to 2300 fleters. 



Ey manipulating the above equation it is possible tc solve 
for tie depth at which the redes should be observed. for 
this particular scerario at a range of 604 meters, codes 
should cccur at integer multiples of 25 meters in depth 
(25n) . Figure 2. 1C shows that these transmission less 
maxima do exist as expected, every 25 meters in depth at the 
stated range. 

Eased on simple physical reasoning one would expect 
refractive lending along the ocean bottom due to density 
differences between tie water and the sediment. Ihe ccctcur 
plot reflects a charge in basic pattern at the interface. 
Ihere is a lending of the TI contours that suggest a refrac- 
tive influence. 

Figure 2.11 shows transmission loss contours from a 
range of 60C meters tc 1350 meters. Again in this region the 
water appears to be dominated by surface reflection. Solving 
for the nodes in this surface interference pattern at a 
range of 1350 meters, it is found that these nodes should 
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appears more accentuated than in the previous figure, lut 
still suggests the influence of refraction at the ccean 
lottcm . 

Figure 2.12 displays transmission loss from 1350 
meters tc past the apex at a range of 2100 meters. In this 
figure, tie dominance of the surface reflection mechanism is 
less obvious and the II patterns become more complicated. It 
is in this region that the influence of trapped normal mode 
propagation can be seen. As discussed earlier, as acoustic 
energy travels up the slope toward the apex, normal modes 
are cut off and energy is transmitted into the bottom. 
According tc adiabatic normal mode theory, medal separation 
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and Zaur 



is rar ce dependent (Graves, Nagel, Uberall, 
and Coppers and Sanders, 1 S80) . The range from the 
wnich the lowest mode is transmitted into the bottom 
calculated using the bellowing equation (Coppens, 
Icanncu, and Kawamura, 1978) : 



, 1 S 75 

apex at 
, can be 
Sanders, 



X = SIN0 c TAN (3 



where: 

X = Dump Distance Of The Lowest Mode 
i) = Naveiergth 
& c = Critical Angle 
£ = Wed ge Angle . 

According tc adiabatic normal mode theory (Kinsler, Trey, 
Coppers, and Sanders, 1982) , in deep water (near the source) 
the lewer normal modes are far above cutoff and the adia- 
batic eigenfunctions consist cf an integer number cf half 
sine laves with zero pressure at both the top and tettom 
surfaces. At the cutoff of each mode, the pressure at the 
bottom must be maximised, resulting in an adiabatic eigen- 
function that contains 1/4# 3/4, and 5/4 wavelengths at the 

respective cutoff distances cf IX, 3X, and 5X for the three 
lowest medes. from figure 2.14 it can be seen that as the 
normal medes travel up the sloping bottom, succesive modes 
are forced into the lettom at distances where a particular 
mode reaches a depth at which it can not longer propagate. 
Also from the figure it is obvious tnat a source set at 
mid-depth can net excite the second mode. Since this 
particular geometry is present in the tank scenario it is 
expected that the energy associated with this second mode 
should ret be seer. Based on this line of physical 

reasoning, modes should be dumped into the sediment at the 
first dump distance, fifth dump distance, ninth dump 

distance and so on. If these modes are dumped as described 
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RANGE FROM APEX 




figure 2.14 Normal Bode Propogation in a Sedge Shaped Ocean. 

then they should he seen in the contour plot as reductions 
cf tr ar emission loss in the bottom at ranges of 1952 meters, 
1760 meters, 1 5 6 € meters, ard so on. 

frcm Figure 2.12 it can be seen that there is an 
cbvicus decrease in tiansmissicn loss at approximately 1950 
meters (dump distance #1) and 1760 meters (dump distarce 
#5). Here is also a less clearly defined reduction in TL 
along tte bottom at 1565 meters (dump distance #9). 

3he final contour plot (Figure 2.13) shovs transmis- 
sion less contours from 1550 meters to 2300 meters. His 
figure is ar extension cf Figure 2.12, intended to emphasize 
how clearly the beam at the first dump distance is defined 
in the plots. The team at the fifth dump distance is also 
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visible hut is not nearly as well defined. From this figure, 
cne can alsc see an indication of a very narrow Lean ir the 
sediment at about I860 meters. Ihis distance corresponds to 
the third dump distance (second normal mode) . Based or. adia- 
batic mode theory the second mode is not expected to he 
excited. However, adiabatic mode theory is only an approxi- 
mation cf normal mode behavior. This approximation of normal 
mode behavior becomes less exact as the bottom slope 
increases and the closer the source is to the apex. Ihe 
appearance cf a narrow beam at the third dump distance indi- 
cates that there is a strong possibility that the second 
mode is present and that the adiabatic approximat icr is ret 
exact with a ten degree hot torn slope. 

Ihe basic features of the IFE contour plots are 
consistent with both physical reasoning and theory. Basic 
surface reflection and bottom refraction occur where 
expected ard behave as anticipated. In the far field, 
trapped normal mode propagation is observed and can he veri- 
fied with simple dump distance calculations. Ihe location of 
beams dumped into the bottom appear consistent with basic 
mode theory. In short, the transmission loss contours indi- 
cate that the IFD is making reasonable predictions cf the 
acoustic field in a shallow water environment. 

5. Ver ifi cation Summary 

It is difficult to say how exact the IFD predictions 
are fer a shallow water environment based on these simple 
verification technigues. As a minimum it can be said that 
the model results are at least consistent with other model 
predictions and expectations based on simple physical 
reasoning. Model estimates are virtually the same as the 
Jenser and Kuperman IF model in regions not influenced ry a 
water/sediment boundary. Close to the boundary the IFD 
results appear to shew greater detail and variation than 
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this IE model. Ihe IID results also appear consistent with 
the general treads la pressure amplitude predictecd by the 
Coppers, Humphries, and Sanders image model. Again differ- 
ences were noted in the small scale structure. Finally, the 
IFD II contours verify well with tasic expectations based on 
physical reasoning and theory. Surface reflection and 
bottom refraction patterns are observed as anticipated. Far 
field normal mode propagation can be verified in the plots 
using simple dump distance calculations. In short, all veri- 
fication methods attempted, fail to uncover any inconsis- 
tency in IFI performance in a shallow water environment. 



4 1 



III. IABOBA ICE I ME ASU R EMENTS 



A. EACKGECUND 

The mayor attempt at appraising IFD performance involved 
comparing model results with laboratory measurements. lie 
shallow water environment modeled in tne tank is very ideal- 
ized; a ten degree sloping sand bottom with an iscvelccity 
water column. Although this scenario appears eiteremiy 
simplistic, it is ore that can he reasonably modeled in tbe 
laboratory and still approximate conditions in an actual 
shallow water ocean environment. The methods used to node! 
and measure the acoustic field are relatively untested. 
Indeed, this attempt at laboratory modeling was performed 
not only to verify HE predictions, but also to see if the 
environment could be successfully modeled in the laboratory. 

E. EXEIFIflINTAl DESIGN 

1 . The Tank 

A fiberglass coated wooden tank was used. The tank 
is 3Gh centimeters ir length, 117 centimeters wide and 55 
centimeters deep. Sand filling the bottom of the tack was 
shaped to form the ten degree sloping bottom, and 
measurements were taken over a range of two meters. Maximum 
water depth in the tank was 35 centimeters. A 100 kHz 
source was placed at mid-channel depth (17.5 centimeters} . 
The layout cf the tark is depicted in Figure 3.1. 
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Figure 3.1 Experimental Tank Set Op 




A slope of ten degrees was selected for several 
reasons. 1c begin with, even though a ten degree slope is 
greater than most ocean bottom slopes, it is still small 
enough tc re considered realistic. Pernaps most important, 
the ten degree slope was selected because given the 
frequency limitations (discussed later in this chapter) and 
range limitations, this wedge angle allowed the source tc be 
placed many (4 1.6) dump distances from the apex. A large 
number ci dump distances was necessary to simulate a distant 
source . 

The bottom material used in the experiment was #30 
fine grace sand. Ihe grain size ranged from 0.15 millimeters 
to 0.70 nilimeteis. Ihe sand was treated with a technique 
used by Eaek (1984) tc remove air from the sediment. This 
technigue used a high speed jet to agitate the sand/water 
mixture tc remove the bubbles and then allowed the sand to 
settle for several days before the experiment was initiated. 

fresh (tap) water served as the medium in the tank, 
lo remove air bubbles, the water was allowed to settle in a 
settling tank before being transferred to the experimental 
tank for use. Ihe water in the tank was periodically treated 
with chlorine bleach to prevent the growth of biological 
material. 

2 . Sig nal Generati ng/ B ece i ving E g u i p ment 

Ihe acoustic signal used for the measurements was 
produced by a function generator, sent through an amplifier 
and then transmitted into the water by a directional trans- 
ducer resonant at 100 kHz. Ihe dimensions of the active face 
of the transducer were 7.0 cm in width and 2.0 cm in height. 
These dimensions resulted in an approximate beamwidth (angle 
from the acoustic axis to the first theoretical null) of 
11.9 degrees in the horizontal and 46.3 degrees in the 
vertical. The narrow horizontal beam minimized reflections 
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from the sidewalls cf the tank, while the wide vertical team 
allowed complete er.soniiicaticn of the channel in the 
vertical dimension. 

lie signal was received ty an LC-5 omni-directional 
hydrophone, sent through an amplifier and filter, and then 
displayed cn an oscilloscope and voltmeter. A schematic 
showing the electronic setup is shown in Figure 3.2. 

A frequency cf 100 kHz was selected for twc prac- 
tical reasons. First, to avoid particle scattering ty the 
sediment, the acoustic wavelength must be at least three 
times larger than the grain size of the sediment (Anderson 
and lietermann, 1966). The largest grain size in the sand 
was C.C7 cm, so that a waveientn of 1.45 cm(100 khz) was 
sufficiently large er.ough to be immune to this effect. 
Second, the 100 kHz frequency and the properties or the sand 
provide a dump distarce that is small enough to allow the 
source tc be positiored many dump distances from the apex. 

Shaping the sand bottom into a ten degree wedge with 
a uniform and smooth interface proved to he a lcng and 
tedious process. ic facilitate this modeling, wooden 
supports (two-by-fours) were mounted along the length of 
both sides c£ the tank. These supports were elevated at one 
end cf the tank to achieve the required ten degree slope. A 
scraping device was constructed with wooden supports and a 
metal scraping blade that extended across the width cf the 
tank. This scraping device consisted of wooden supports 
along the top that reached across the tank and could be 
pulled along the elevated wooden supports on both side of 
the tank. This scraping device was pulled along the supports 
repeatedly until a smooth slope of ten degrees was sculp- 
tured frcm the sand. 

Fcles had tc be drilled into the metal scraping 
Hade because when a solid blade was used in the shallow 
portion cf the slope, water trapped behind the blade was 
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1. ) LC-5 HYDROPHONE. 

2. ) AMPLIFIER (HP 465A). 

3. ) ELECTRONIC FILTER (SK 302). 

4. ) OSCILLOSCOPE. 

5. ) VOLTMETER (HP 400D). 

6. ) FUNCTION GENERATOR WAVETEK 116. 

7. ) FREQUENCY OSCILLATOR (GR 1310). 

8. ) AMPLIFIER (HP 467A). 

9. ) FREQUENCY COUNTER (HP 5233L). 

10.) 100 KHZ TRANSDUCER. 



Electronic Equipment Schematic. 
4 6 



figure 3.2 



forced beneath the scraper gouging the smooth hotter. This 
hotter modeling techrigue was slow because, once a pass was 
made ever the better, the water became turbid and it was 
then impossible to see the bottom. Wnen the bottom was not 
visible, it was impossible tc see where further smoothing 
was necessary until the water settled several hours later. 
In addition, as this smoothing process continued a silty 
residue became separated from the sand and settled cut on 
top of the sand. Ibis residue would be easily resuspended 
and eventually had tc be removed using a water syphon. 



C. MfASGEEbENT EEOCIIUEES 

Measurement of the pressure field within the water was 
done by lowering the receiver in depth at specific ranges of 
interest. Ihe receiving hydrophone was attached to a pair of 
micrometers at right angles to one another, that was in turn 
bolted tc a board which spanned the width of the tark. Once 
the beard was placed close to a range of interest, cne 
micrometer was used to give fine adjustments in range and 
the ctker in depth. 

Ihe measurements were subject to both an accuracy and a 
precision error. On a given day, with the water level fixed 
and tbe cross-tank support set at a particular place in the 
tank, it was possible to position the receiving hydrophone 
with an accuracy in depth and range of plus or minus 0.06 
centimeters {one turr of the micrometer) . lo prevent the 
sand 'inland* of the apex from drying out when not taking 
measurements, enough water was added to the tank after each 
data run tc keep the sand completely submerged. Ihe next 
time measurements were taken, water had to to be removed 
from the tank tc reestablish the beach. Because of these 
small changes in the water level the horizontal position of 
the bead was subject tc a precison error estimated to be 
withir plus or minus cne centimeter. 
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Ike final decision associated with the a easure merts 
centered cn whether tc use a triggered pulse cr a continuous 
wave(CK) signal. Kith a triggered pulse it was possible to 
distinguish the received signal iron interference caused by 
■teflecticns off the side of the tank. On the other hare, by 
using a triggered pulse there was a possibility that tne 
pulse length was not long encugh for the acoustic energy 
associated with paths reflecting off the top and bottom of 
the water column to overlap the direct path from the scarce 
to receiver. A CW signal would avoid potential pulse length 
problems, but it would be impossible to distinguish between 
the actual signal and interference. The use of a source with 
a narrow horizontal team reduced the effects of reflections 
from the side walls, tut there was still the possibility of 
interference from side lobes reflected from the sides. 

It was necessary to determine the best means tc take 
measurements. This was done by taking measurements at the 
same location with different pulse lengths to determine the 
required pulse length to give consistent results and then 
comparing these results to those obtained with a CW signal. 
Tne third dump distance (14.4 cm from the apex) and just 
past the tenth dump distance (50.0 cm from the apex) were 
chosen, and measurements were taken for pulsed signals of 64 
and 256 cycles, and a CW signal. These measurements can be 
seen in figure 3.3 through figure 3.5. In the figure 
depicting results at the third dump distance , the CW 
measurements are shown as the solid curve, the 64 trigger 
cycle results are displayed as circular points, and the 256 
trigger cycle results are depicted as triangular points, for 
clarity, the measurements taken just past the tenth dump 
dist a rce axe shown ir two figures. Figure 3.4 compares CW 
results (solid curve) with the 64 trigger cycle signal 
(circular points). figure 3.5 compares CW measurements 
(solid curve) with the 256 cycle results (circular points) . 
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figure 3.3 lulse length Analysis at 3. OX. 
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figure 3.4 lulse length Analysis at 10. 4X. 
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figure 3.5 Eulse length Analysis at 10^4X. 
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Ihe pressure amplitude values were normalized by dividing 
each pressure value by the maximum pressure value ir the 
field. Ihe depth values were also normalized by dividing 
each depth ty the maximum depth. 

Ihe results in figure 3.3 show good agreement between 
the three curves. In particular, the results for the CW 

signal ard the 256 cycle signal are very similar. Ihe 64 

cycle signal shows the same behavior as the CW sigral, but 
differs in magnitude below mid-depth. 

The results at 5C.0 cm frea the apex appear much more 
complicated than at the third dump distance. Figure 3.4 
compares results acheived with the Cw signal ana the 64 
cycle signal. From the figure it can be seen that the two 
sets cf results shew good agreement. The measurements 
depicted in Figure 3.5, which compares the 256 cycle signal 
and the CW signal show poorer agreement, although the 
general shapes cf both curves remain similar. 

Ihese results indicate that interference from the side 
walls is sufficiently small sc that it is possible tc make 
measurements with a Cf signal to at least ten dump distances 
from the apex. Long pulses could also be used, tut the 
difficulty cf making voltage measurements with an oscillo- 
scope compared to reading a voltmeter, dictates that 
measurements should he made with a CW signal. 
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IV. COD II EESUII CCCIAEISONS SITH 1ABOBATOE Y aE^SUSEMEEIS 



A. IB1ECECCTION 

I ED model predictions and laboratory measur erects were 
cbtaired fcr comparison in both a general and detailed anal- 
ysis. The general analysis compared results every rive 
centimeters from the beach to a range of 50.0 cm. Ires the 
leach. The measurements were spaced to give results at 
a r prcximately each cf the first ten dump distances. The 
detailed analysis compared results every centimeter frcn 3.0 
to 11.0 cm from the apex. These measurements were taken to 
observe the sensitivity of the pressure field tc small 
changes ir range. 



E. TEE GENEEAL ANALYSIS 

The general anlysis compared IFD results and laboratory 
measurements at approximately eacn of the first ten dump 
distances (every five centimeters from the beach). These 
measurements were obtained by fixing the receiving hydro- 
phone with respect tc the cress board and then moving the 
hoard out in increments of five centimeters to measure the 
desired field as a function cf depth at each range. The 
received CK signal was read on the voltmeter. 

The comparisons are shown in Figure 4.1 through Figure 
4.10. In the figures, the IFD predictions are displayed as a 
solid line while the experimental measurements are Shewn as 
circles connected by a dashed line. Each pressure value was 
normalized to unity hy dividing it by the maximum pressure 
value fcr the respective curve. The depth values were also 
normalized ly dividing each depth by the maximum depth in 
the water column. 
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figure 4.1 Comparison cf Besults at 1.0X. 
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figure 4 .<2 Comparison or Besults at 2- IX. 
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figure 4.3 



Comparison of Results at 3. IX 
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figure 4.4 Ccaparisou o£ Results at 4.2X. 
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Figure 4.5 Comparison of Eesults at 5.2X. 
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figure 4.6 Comparison of Results at 6.2 X. 
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figure 4.7 Comparison of Eesults at 7.3X. 



60 



H/Z 



NORMRLIZlD RMPLITUDl ( F/PMRX ) 

0.0 0.2 0.4 0.6 O.B 1.0 



1 .2 




figure 4.8 Ccmparisoc of Results at 8.3X. 
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Comparison 
6 2 



figure 4.9 



of Besults at 9.4X 




Figure 4. 10 



Comparison 



of Eesuits 



at 



10. 4X. 
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Ir general, the pressure patterns predicted by the JED 
and these measured, both become more complicated as the 
range frem the beach increases. At all ranges, there is 

qualitative agreement in the scale of the predicted features 
and the scale of the measured features. Quantitative agree- 
ment is lacking. This agreement in scale but not detail, 
suggests that the phases ameng the normal modes predicted by 
the mcdel do not accurately reflect the experimental 
situa tier. 

C. Til LET AILED ANA1ISIS 

The detailed analysis compared IFD values witn experi- 
mental measurements every centimeter from 3.0 cm to 11.0 cm 
from the teach (0.7X tc 2. 3X) . The laboratory measurements 
here taken by fixing the hoard at one location and then 
using the micrometer to adjust the receiver to the desired 
range. The received CJi signal was read from the voltmeter. 

These results are depicted in Figures 4.11 tnrough 4.19. 
In the figures, the IFD predictions are shewn as a solid 
line and the experimental measurements are displayed as 
circles. Each pressure value was normalized to unity by 
dividing it by the maximum pressure for the respective 
curve. The depth values were normalized by dividing each 
depth hy the maximum depth at the particular range. 

There is qualitative agreement in the scale of the basic 
features fer all ranees, but quantitative agreement is ret 
observed. The IFD patterns change very little throughout the 
analysis, while the measured values change more rapidly 
(especially past 2. OX). The results at 2. IX and 2.3X 
(Figure 4.16 and Figure 4.19) snow that the pressure field 
can change fairly significantly over a range as short as cne 
centimeter. At those ranges where tnere is poor agreement 
between results, there is an indication of phase 
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figure 4.11 Comparison of Results at 0.7X. 
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Figure 4.12 Comparison of Results at 0. 8X. 




figure 4.13 Comparison of Results at 1. OX. 
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figure 4.14 Comparison of Hesults at 1.3X. 
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NORMRLIZlD RMPLITUDE ( P/Pr.RX ) 
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figure 4.15 Comparison of Results at 1. 5X. 
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figure 4.16 Coapariscn of Results at 1. 7X. 
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NORMALIZED AMPLITUDE! P/PMAX ) 
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figure 4.17 Comparison of Results at 1. 9X. 
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figure 4-18 Comparison of Results at 2. IX. 
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figure 4-19 Comparison of Besults at 2- 3X. 
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interference (Figure 4.18) . For distances less tha 
theory predicts that for a source at mid-depth, the 
propagating mode should experience interference from 
the evanescent tails of higher modes. Therefore, 
region of the detailed analysis, one expects 
interference effects. The experimental results in 
region however, shew rather significant interference 
interference suggests unsuspected propagating mode 
present, or that tie evanescent tails are larger 
expected. From the trend in the curves it appears as 
phase interference is not a ractor from 1.5X inward 
the teach (Figures 4.11 through 4.15) . 
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V. COftClOSIONS/EECOMMBNDATIONS 



A. CCNCICSICNS 

1 . I erfor ma ace cf the I f D M ode l 

Ihis analysis of the IfD acoustic model did not 
uncover any major failures cx model performance in a simpli- 
fied shallcw water environment. But this is not to say that 
the model consistently and accurately performs in such an 
ccean scenario. Although comparison of IFD predictions with 
two cttei mcdels and with simple physical, reasoning did not 
uncover any inccnsist encie s in performance, the agreement 
between IIE values and laboratory measurements is insuffi- 
cient to give complete confidence in the performance cf 
either tie model or the experiment. 

Ihere is reasonable agreement between the scale of 
the features predicted by the IfD model and the experiment. 
Eoth the model and the measurements show increasing 
complexity as range is increased from the beach. But despite 
the similarities , there is poor guantitative agreement in 
results. Cre possible cause for the differences may be that 
the phases of the normal modes predicted by the model are 
extremely sensitive to minor irregularities in the shape of 
the interface and the acoustic properties of the bottom, 
freliminary work (LeSesne, 198h ) , suggests that the phase 
relationships between modes is strongly dependent on the 
distance of the source from the apex even at great ranges. 
Conse cuertly , it appears that the collective influence cf 
the normal modes is dependent upon careful geometric control 
cf the experiment. 

Ihe detailed analysis arcund 2X, reveals that large 
phase interferences occur where only one propagating mode is 
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expected. Inis indicates the existence of an extended 
evanescent tail within cutoff cf the higher modes. Kavemcra 
and Icanncu (1 978), noted that the apparent phases cf the 
evanescent tails are extremely sensitive to the details of 
the ervixcr. nen t. The results indicate that the tail dees not 
decay quickly, and its influence is pervasive (Figures 4.2 
and 4.2) . From the trend in the curves, it appears this 
phase interference is not a factor any closer tc the teach 
than 1 . 5 X (figures 4.11 through 4.15). 

2 • Kc d eli n g /Measur eme nt Frcc edu res 



The lab oratory measurements represent an initial 
attempt at modeling a:, idealized shallow water environment. 
The experimental techniques are not without problems. All 
sets cf measurements were repeated and snowed good agree- 
ment. The reproducibility cf the measurements suggest that 
random errors have been minimized. However, systematic 
errors remain which also contribute to the dr sc repancies 
between predicted and measured values. Some of these errors 
derive from equipment, such as the fact that the size cf the 
hydrophone is of the same order or magnitude as the scale cf 
variatiers in the pressure field. In addition, the environ- 
ment in the tank may net be sufficiently close to the ideal 
environment assumed by the model. 



E. EECCEEEEDATIGNS 

fritter study and verification of tne IFD computer model 
is recommended. A comparison cf IFD predictions with ether 
models net analyzed in this study, may uncover the cause of 
the inconsistency when compared to measured values. A 
detailed study of the phase interaction of the normal modes, 
although extremely complex, may also offer insight into the 
If D performance. 
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It appears that the modeling of an idealized shallow 
hater environment in the laboratory is of vaiu^ xn acoustic 
model verification. tespite the experimental difficulties, 
there has qualitative agreement in the basic features 
letween ircdel predictions and laboratory measurements. Eut 
the experiment represents only an initial probe into the 
modeling/measurement techniques of a shallow water environ- 
ment. further refinement of these techniques and the experi- 
mental equipment may result in tetter -uantitarive 

between model predictions ai .. _ . rur at _ 

.. casat ments. 
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APPENDIX A 



BEVISED IFD PECGEAtf LISTING 
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APfENIIX C 

BURNING THE CCNTOUB EICT ON THE NPS COflPXTEB 
A. INTECDUCTION 

This Appendix describes a procedure for running the Ti 
contcur plot on the NES computer. Detailed instructions tor 
running the IFD program can he found in Jaeger {1 9 83} . 

E. CCEY1NG THE FIXE FOE USE 

Cnee the TI values have been generated by the IFD 
program ail that is needed to produce a II contour plot is 
the file EICTS FCE1BAN. This file is shewn in Appendix E and 
can he copied from a computer account maintained hy the 
Underwater Acoustics Curriculum. To link with this account 
and chtain a copy of the file, the user should proceed as 
fellows: 

(1) leg on terminal. 

(2) Enter: CP LINK 0160P 1S1 195 ER . 

(3) Wten prompted for the read password enter: UX . 

(h) Erter: ACC 195 C . 

(5) Enter: COEY PICTS FOE1EAN C = = A . 

At this point the PICTS FOEIEAN file should reside on the 
user’s A disk. 

C. BURNING THE EBOGEAN 

Eefcre running the program the user must obtain the TI 
values from the IFD program. These TL values must he sent to 
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CSC 



a data disk in cider ior the fICIS program to he able tc 
them. This can he dcre hy placing a WRITE statement in the 
IFD program that sends the values to a data disk. The modi- 
fied HE depicted in Appendix A uses this technique and can 
he used as a guide. 

Cnee tie data disk has teen created, the user mist 
assign temporary disk space (1DISK) to give the program 
sufficient room to generate the 11 contours. For mere infor- 
mation cn how to assign temporary disk space, see iPS 
Technical Note TN-VH-C1 which is availitle in the computer 
consultant’s office. 

With these initial steps completed, all that remains is 
tc compile and run tie program. The program can he compiled 

ty: 

Enter: FCETGI PLOTS . 

The program must he run at the TEK618 graphics terminal 
under DISSPIA. This can be dene hy: 

Enter : DISSPLA . 

The user will then he prompted for the compiled Fortran 
program name and the file definitions for the data disk, 
hefore the program will run. 
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APPENDIX D 

SOUfiCE IFPTH SENSITIVITY ANALYSIS 



During the laboratory experiment, two supplemental sets 
cf measurements were taken to obtain an indication of the 
sensitivity of the pressure amplitude to changes in the 
source depth. The first set of measurements was obtained by 
varying the source depth with the receiver fixed at 1.0X 
<4. 8 cm from the beach) , and lowered in depth to the bottom, 
for tie second set of measurements. The receiver was fixed 
at the third dump distance (14.4 cm from the beach), and 
measurements were taken with the source fixed at 5, 7, 9, 
11, 11, 15, 17, 19, 2 1, 23, 25, 27, 29, and and 31 cm from 
the surface. 

The results of the first analysis are shown in Figure 
C. 1. In the figure, pressure was normalized by dividing by 
the maximum pressure and depth was normalized by dividing by 



Although theory predicts only 
the results show 



the depth cf the water column 
one propagating mode at this distance, 
modal phase interference. 

The results cf the second analysis are shown in Figure 
D.2 through Figure D.6. In the figures, each curve repre- 
sents a set of measurements taken with the source fixed at a 
specific depth. Fox convenience, more than one curve is 
shown on a given figure. All pressure values were normalized 
by dividing by the the maximum pressure in the field. The 
depth values were normalized by dividing by the maximum 
water depth. 

Ir general, at 2. OX there appears to be modal phase 
interference which is influenced by source depth. Giver the 
roughness of this analysis, the details of this interference 
are obscure. However, the modal interference does net appear 
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to he incor. sistect kith either the placement of the source 
cr with the previous experi nectal measurements. 
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figure D« 1 Source Sensitivity Analysis at 1.0X. 
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figure D.2 Mearur eaents with Source Depth of 5, 1, and S Cm. 
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liguxe D.3 fleasureieuts with Source Depth of 11 # 13, and 15 
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Measureaents with Source Depth of 17, 19, and 21 
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Figure D.4 
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Eiguxe D.5 Heasureaents with Source Depth of 23, 25, and 27 Co. 
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figure D.6 Bea£urei€cts with Source Depth of 29 and 31 Cm. 
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